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Abstract 

Even though the Standard Model with a Higgs mass m H = 125 GeV possesses no bulk 
phase transition, its thermodynamics still experiences a “soft point” at temperatures around 
T = 160 GeV, with a deviation from ideal gas thermodynamics. Such a deviation may have 
an effect on precision computations of weakly interacting dark matter relic abundances if 
their mass is in the few TeV range, or on leptogenesis scenarios operating in this temperature 
range. By making use of results from lattice simulations based on a dimensionally reduced 
effective field theory, we estimate the relevant thermodynamic functions across the crossover. 
The results are tabulated in a numerical form permitting for their insertion as a background 
equation of state into cosmological particle production/decoupling codes. We find that Higgs 
dynamics induces a non-trivial “structure” visible e.g. in the heat capacity, but that in general 
the largest radiative corrections originate from QCD effects, reducing the energy density by 
a couple of percent from the free value even at T > 160 GeV. 
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1. Introduction 


The excellent performance of the Standard Model (SM) in describing LHC data suggests that 
the SM represents a precise description of nature up to energy scales of several hundred GeV. 
If this observation continues to be confirmed by future runs, one consequence is that the Higgs 
phenomenon set in smoothly, i.e. without a phase transition, as the Universe cooled down to 
temperatures below 160 GeV [1]—[4] (cf. ref. [5] for a review of refined investigations). This 
would mean that the bulk motion of the plasma did not depart from thermal equilibrium and 
therefore did not generate cosmological relics. 

Nevertheless, the thermodynamics of the SM could still have left a “background” imprint 
on other non-equilibrium physics that just happened to be going on. One example is that the 
B+L violating rate switched off rapidly around the crossover [6], and therefore determined 
which fraction of a given lepton number being produced around these temperatures could be 
converted into baryons (cf. e.g. ref. [7]). Another is that if Dark Matter particles decoupled at 
around this period, even small features in the equation of state could have had an impact [8]. 
(Similar considerations can be carried out for the QCD crossover, cf. refs. [9, 10] and references 
therein.) A rough estimate for the decoupling temperature of weakly interacting particles of 
mass M is T ~ M/25, so that we may expect the largest sensitivity for M ~ 4 TeV or so. 

The purpose of our study is to estimate the equation of state of the SM around the elec- 
troweak crossover. We do this through perturbative computations extending up to 3-loop 
level, as well as by making use of existing lattice simulations within a dimensionally re¬ 
duced effective theory. Only the temperature is treated as a non-trivial canonical variable, 
with chemical potentials associated with conserved charges such as the baryon minus lepton 
number or the hypercharge magnetic flux set to zero. 

The paper is organized as follows. After defining the basic observables in sec. 2, we derive 
a “master equation” in sec. 3 which relates the trace of the expectation value of the energy- 
momentum tensor to three ingredients. The first ingredient, scale violations through quantum 
corrections, is addressed in sec. 4. The second ingredient, the expectation value of the Higgs 
condensate, needs to be evaluated up to the non-perturbative level, and this is achieved in 
sec. 5. The third ingredient is related to vacuum renormalization and is discussed in sec. 6. 
All ingredients are put together in sec. 7, where we also present phenomenological results. 
Section 8 contains some conclusions as well as a discussion of the theoretical uncertainties of 
the current analysis. Readers not interested in technical details may start from sec. 7. 

2. Basic setup 

We consider the Standard Model with a Higgs potential of the form 

5C e = + A B (0 t 0) 2 , (2.1) 
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where v‘ B , A B are bare parameters (the corresponding renormalized parameters are denoted by 
v 2 , A), and C E is a Euclidean (imaginary-time) Lagrangian. Denoting by g 2 € {A, h 2 , g \, g 2 , g$} 
a generic coupling constant 1 , we concentrate on the parametric temperature range 


> d 

r^j o 

9 2 


( 2 . 2 ) 


in the following. At the lower edge of this range, the effective Higgs mass parameter of the 
dimensionally reduced theory [11, 12], m 2 [13], is assumed to satisfy 

\m 2 \ ~ \-v 2 + g 2 T 2 \ < (2.3) 

7r 

but it can have either sign. If it is negative, we may expect to find ourselves in a “Higgs 
phase”. Within a gauge-fixed perturbative treatment, this would imply that the Higgs field 
has an expectation value v 2 ~ —fhj/A > 0. Within the range of eq. (2.3), only relatively 
small expectation values v<g^T can be considered, where we counted A ~ g 2 . 

It is well known that when momenta in the range | m 2 \ ~ ( g 2 T/it) 2 are considered, which is 
a special case of eq. (2.3), then the dynamics of the system is non-perturbative [14, 15]. There¬ 
fore the dynamics needs to be treated with lattice methods. However, non-perturbativity is 
associated with particular modes only, and can be captured with a dimensionally reduced 
effective description [11, 12], Even though the construction of this theory is perturbative, we 
expect to obtain a description accurate on the percent level within a weakly coupled theory 
such as the SM. This estimate is based on analyzing the influence from higher-order operators 
that are truncated from the dimensionally reduced theory (cf. sec. 5.4 of ref. [13]) and, more 
concretely, from a comparison of non-perturbative results for the location of the endpoint in 
the phase diagram of an SM-like theory based on the dimensionally reduced description and 
on direct 4-dimensional lattice simulations [16, 17]. (It would be interesting to repeat the 
latter type of a comparison for the physical value of the Higgs mass.) 

The basic observable that we consider is the thermodynamic pressure of the SM. The 
pressure can formally be defined through the grand canonical partition function Z as 


Z = exp 


(M^). 


(2.4) 


where the thermodynamic limit V —>■ oo is implied, and p B denotes the bare result. The 
computation of p B has previously been considered in ref. [18] for the case m 2 ~ g 2 T 2 , and in 
ref. [19] for 777.3 ~ g 3 T 2 /ir. Following refs. [18, 19] 2 (which were inspired by ref. [20]) we write 


Pb( T ) = P e { T ) +Pm( T ) +Pg( T ) > ( 2 - 5 ) 

1 Here ht denotes the top Yukawa coupling and pi,i? 2,<?3 are the gauge couplings related to the Standard 
Model gauge groups U Y (1), SU L (2) and SU C (3). 

2 We agree with most of the results in these papers, however disagree with the renormalization condition 
for p and with certain technical details, cf. appendix B. 
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where p E collects contributions from the momentum scale k ~ n tT, p M those from k gT , 
p G those from k ~ g 2 T /7r. Rephrasing terminology often used in the QCD context, we refer 
to the effective theory contributing to p M as ESM (“Electrostatic Standard Model”) and to 
that contributing to p G as MSM (“Magnetostatic Standard Model”). 

As defined by eq. (2.4) the pressure is ultraviolet divergent. We renormalize it by assuming 
that the pressure (and energy density) vanish at T = 0. The renormalized pressure can then 
be written as 

p(T)=p B (T)-p B ( 0) . (2.6) 

Our goal is to determine the dimensionless functions p(T)/T 4 and T<9^{p(T)/T 4 } up to 
0(g 5 ). For the latter quantity, it turns out that the contribution of the softest momentum 
modes needs to be treated non-perturbatively in order to reach this precision. 


3. Master equation 


Consider a normalized form of eq. (2.6), p/T 4 = \p B (T) — p 0B ]/T\ where p 0B = p B (0) de¬ 
notes the bare zero-temperature pressure. Unless stated otherwise ultraviolet divergences are 
treated through dimensional regularization. In the difference of eq. (2.6) all 1/e poles cancel, 
so that the bare expressions can be replaced with renormalized expressions in the MS scheme 
0 Pb Pr, Pob -> Por) : 


p{T) = p R (T;p,is 2 (p),g 2 (p))-p 0R (p,is 2 (g),g 2 (p)) 

'J'i 


(3.1) 


Here p denotes the MS renormalization scale. 

We envisage that at some temperature (T = T 0 ) on either side of the crossover, say 
T 0 <C 160 GeV or T 0 S> 160 GeV so that \rh 2 \> g 3 T 2 /tt, p/T 4 can be determined by a 
direct perturbative computation. The task then is to integrate p/T 4 across the electroweak 
crossover: 


_) P(T 0 ) _ [ T i dT d fp(T)\ 

Tf T 4 J Tq T dT \ T 4 J ‘ 

So, we need to compute the logarithmic temperature derivative of the dimensionless ratio 
p/T 4 . The result is non-zero because of the breaking of scale invariance, either by the 
explicit mass term is 2 , or by quantum corrections related to running couplings. 

Making use of standard thermodynamic relations, we note that 



T d f p(T) ) 
dT\ T 4 J 


e(T)-3p(T) 

J'A 


= A(T) , 


(3.3) 


where e denotes the energy density. In the context of QCD this quantity is referred to as the 
trace anomaly, and has been studied with lattice methods. Were it not that chiral fermions 
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(in particular the top quark) are essential for the physics considered, the problem could in 
principle be studied with full 4-dinrensional lattice simulations here as well. In the present 
paper, we circumvent the problem of chiral fermions by using lattice input only for the Bose- 
enhanced infrared degrees of freedom, treating fermions within a (resummed) weak-coupling 
expansion. 

Now, because of dimensional reasons, we may rewrite eq. (3.1) as 


P(T) _ 


P V (fi) 2 lr . 


~ J%(p,V 2 (p),g 2 (p)) • 


We note from the Euclidean path integral representation, viz. 


(3.4) 


2 = 


exp{-^d 3 - 2 e x^ VT dr£ E } 


that the bare pressure obeys (here p B = p B /T 4 


(3.5) 


dp E 


[Zm(4>* 4>)]l 


(3.6) 


d{v 2 {p)/T 2 ) T 2 

where we wrote u 2 = Z m u 2 (p). 3 It is now straightforward to obtain the following relation: 


,_d_ f P(T) \ = dp B 


dT( T 4 J 9 In {fi/T) 


2v 2 (p) [Z m (4^<j>)\ R + 4p 0R 


rp4 




(3.7) 


Here eq. (3.6) has been replaced by its renormalized version. 

It can be observed from eq. (3.7) that three ingredients are needed: the determination of 
explicit logarithms appearing in p R (“breaking of scale invariance by quantum corrections”); 
the temperature evolution of the Higgs condensate (“explicit breaking of scale invariance”); 
as well as the vacuum term needed for renormalization. We discuss the first of these in the 
next section, the condensate in sec. 5, and the vacuum term in sec. 6. The results are collected 
together and evaluated numerically in sec. 7. 


4. Effects from scale violation by quantum corrections 

The first ingredient needed in eq. (3.7) are logarithms of p/T that are induced by loop 

dv 

corrections. There are two ways to extract g i n (g/r) f rom re f- [18]: either by reading logarithms 
from the explicit expressions for the various coefficients given there, or by deducing them from 
the running of the couplings. Making use of the notation in eq. (2.5), with p = p/T 4 , the 

3 The renormalization factor reads Z m = 1 + [2A + h% — | (gf + 3^2)] + CKff 4 ) but is not separately 

needed, because it always appears in the combination Z m ((j>'(j>). 
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part p G contributes at 0(g G ). From p E + p M we need terms up to and including (D(g 5 ). The 
expansion reads 


Pe “I" Pm a El T 9l a EB “I" 92 a BA T 93 a EC T '^EA T h-f-CX EY + r ^2° L Ev 


(4.1) 


+ 


Yi + 0 JU 1 


( m ^) 2 + 3 ( toe 2 (0) ) 2 + 8 (? nE 3 (0) ) 2 


+ ... , 


( 47 t) 2 T 4 \e ' ~ El/l/ J ' 12itT 3 
where the leading-order Debye masses read (n G = 3 denotes the number of generations) 


m 


2 ( 0 ) _ 


“ ( 6 + 


5 n c 


)s\T 2 . 


m 


2 ( 0 ) _ 


(jU^) 5 2 T 2 , mii 0} = (l + ^) 5 2 T 2 . (4.2) 


All coefficients in eq. (4.1) apart from a Evv are scale independent; its value differs from that 
in ref. [18] because of our different renormalization condition: 


a 


E vv 


d F n s 




(4.3) 


where d F ,n s are specified in eq. (B.l). Given that the pressure as a whole is scale independent, 
running from the couplings must cancel against explicit logarithms; therefore, 


&\Pe 4~ Pm]r 
dln(p/T) 


_ d g\ _ dg 2 _ dg 2 

" EB n ' d/T “ " BA d^ “ " EC A d^ “ " EA 


_dA _ dtif a El/ _di /2 


+ 


da E 


(47r) 2 T 4 ^ dp 


E 

1=1 


djVn F } 

8vrT 3 


!>- 


. d ml (0) 
dp, 


+ 0(g e 


(4.4) 


where d\ = 1, d 2 = 3, and d 3 = 8 . The runnings can be read from eqs. (7)-(12) of ref. [18]. 
Putting together, we obtain 


d[j>E 4~ PmIr 
9 In (p/T) 

A 1 (T) 


A AT) 

1 


(4.5) 


198 + 141n G — 20n G 4 | 266 + 163n G — 40n 2 4 
93 d" r>oo 92 


(47r ) 2 54 

144 + 375n G + 1000n] 


7776 


O 4 _9bl_ h 2(_ 
yi 32 * V 32 


288 
7 hf 


— A(^A 

1 

(47r) 3 

- 4g 3 


+ -A - 


9i +3 ff 2 \ n 


+ 


y 2 


(h f 2 + 2A - 


M 

6 

9i + ‘igl 


329|(l + y) § (h-^) +1 2 9 | 


5 Uq 

6 + 3 


15 gl 85 g\ 
64 576 

2i/ 4 
J'i 

I /43 
V24 


% 

3 


1 5n G 
6 + 9 


J_ 1 5n Q 
24 + 9 


0( 5 6 


(4.6) 


This expression is renormalization group (RG) invariant up to 0(g G ). The numerically largest 
corrections originate from terms involving the strong gauge coupling g 2 . 
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5. Effects from the Higgs condensate 
5.1. Outline 


The Higgs condensate can be obtained from eq. (3.6). For a non-perturbative evaluation, we 
make use of the lattice simulations in refs. [21, 6 ]. This means that we split the pressure into 
contributions from various momentum scales like in eq. (2.5): 


*•<*> = £ + & 


M _|_ &Pq 


du 2 


(5.1) 


Given that Z m (<jr(j)) is multiplied by u 2 in eq. (3.7) and that we assume v 2 ~ g 2 T 2 , we only 
need to determine Z m (<$(\)) to 0{g 3 ). That said, it turns out that computations going beyond 
those in refs. [18, 19] are needed. 


5.2. Perturbative contributions 

The first term of eq. (5.1), the contribution from the “hard” scales k ~ 7 tT, can be directly 
extracted from the results of ref. [18]: 


dp E 

du 2 


— T < ay, + 


d~ ff2 ^EA v + A «EA, + h: t a E 

(47r) 2 


+ 


2 v 2 
(47r) 


- 


+ 0{g i ). (5.2) 


Inserting coefficients from ref. [18] and from eq. (4.3), this contains 1/e-divergences in the 
terms proportional to g \, g\ and u 2 : 


a Jt = 


1 


du 2 


6 \ (4tt ) 2 [2 


x (gf + % 2 ) ( - + 3 In -j— + 7 e + - + 


5 , 2C'(-1) 


+ 6 h t In 


87 tT 


+ 




1 


+ 7e + 


P 


e 4irT 

12\( In 


C(-i) 


v 47 tT 


+ 7e 


(4tt ) 2 \2e 4 t tT 


+ ln-^+ 7E +0(5 4 )- 


(5.3) 


The second term of eq. (5.1), the contribution from the “soft” scales k ~ gT, originates at 
0(g 3 ) from the dependence of the effective mass parameters on u 2 \ 


dp M 

du 2 


dmli dp M dm | dp^ 6 

du 2 dm.du 2 dm 2 


(5.4) 


The first two terms give contributions that can be extracted from [19]: the Debye mass 
parameters depend on u 2 as 


2 2 

m E1 = g 1 


T 2 (3' E1 - 


(4?r) 


• PL, 


2 2 

) m E2 — 9 2 


T 2 /3 E1 - 


(4tt) 


r/3 E 


+ ... , 


(5.5) 
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d Pm 

d m 2 


jcQ? + HD'H 


Figure 1: Processes representing dp M /dm §. The filled blob is m§ </>'</>; dashed lines are scalar prop¬ 
agators; solid lines are adjoint scalar propagators; and wiggly lines are gauge fields. All 1-loop and 
2-loop diagrams and all other 3-loop diagrams vanish in dimensional regularization. 


where /3 E1 and j3 E1 are as in eq. (4.2). Differentiating the term on the second line of eq. (4.1) 
and inserting the values of the coefficients from ref. [19] yields 


\ - dml, 

^ dv 2 
1=1 


d P M 

dmli 


T 2 

(47t) 3 


q a I 5 . n G 
92 ( 6 + 3 





+ 0(5 4 )- 


(5.6) 


The last term of eq. (5.4) also contributes, but this contribution cannot be extracted from 
ref. [19]. The reason is that if we consider p M without a derivative, then contributions involving 
m 3 are of 0(<? 6 ) for Imjl <g 3 T 2 /-K. However, the derivative of m 2 is larger than m 2 itself: 
v 2 dm 2 /dv 2 ~ v 2 ~ g 2 T 2 . Therefore terms that are of higher order in p M do contribute to 
the trace anomaly. 

The computation of dp M /drn 2 is simplified by the fact that since the softest physics is 
captured by a study of p G , the Higgs and gauge fields can be treated as massless in the 
computation of the matching coefficient p M . Then most diagrams, in particular all diagrams 
which do not contain at least one adjoint scalar propagator, vanish in dimensional regular¬ 
ization. The leading non-zero contributions, which are the only terms needed at 0(g 3 ), are 
those given in fig. 1. Reducing the diagrams to master integrals and evaluating them with 
standard techniques (cf. e.g. ref. [ 22 ]), we obtain 


dm, 2 dp M 
dv 2 dm 2 


Tl [jL('_?_ 12 hl -g- + g') + _^. + i ' 

16(4vr) d L^b 2 V e 2m E 2 3/ m E1 m E1 + m E2 _ 


(5.7) 


5.3. Non-perturbative contribution 


The third term of eq. (5.1), the contribution from the “ultrasoft” scales k ~ g 2 T/-K, can be 
expressed as 


dp G 

dv 2 


(9m| dp G 
dv 2 dfh\ 


+ 0(g 8 ), 


(5.8) 


where m 2 is the Higgs mass parameter within MSM and the error comes from partial deriva¬ 
tives with respect to the other effective couplings of MSM. The dependence of mi on v 2 is 
known up to 0(g 2 ) [13], with the leading term reading dm^/dv 2 = —1. The condensate can 
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be expressed in dimensional regularization as 


dm\ 


(Ami + 3gl 
(47t) 2 


~(r + ln_ i~ 

V 4e 5 m2 


+ {4^ 4) 3d id] 


2 ) 

M2 / 


(5.9) 


where the argument of (4>^4 >} 3c j indicates the renormalization scale used within MSM. As 
discussed below eq. (2.3) the condensate is parametrically <gT 2 in our power counting, and 
therefore we keep the correction of 0(g 2 ) in the coefficient of (^(/>) 3d , even though we do not 
keep it in the ultraviolet part where it is unambiguously of 0{g !i ). For future reference we 
express the ultraviolet part in terms of the couplings of the full theory, by inserting 4 


9m2 

2 

9 Ml 


glT 


9%T 
247t m E2 


glT + 0(g 4 ). 


^1 + 2e In 



+ C ) (5 4 ) > 


(5.12) 

(5.13) 


Thereby the ultrasoft contribution from k ~ g 2 T/t: becomes 


dm 2 dp G 
du 2 dm 3 


= u + 
+ 


2(47r ) 2 
(gj + Sg^T 2 


(. 9i + 3 ff2 - 8A) ln (^ T ) - 4h t ln (^r) } (<^>3 d(9: 


(47r) 2 
9 4 2 T 3 


1 ,i 9 
— +ln — 

4e 5 m2 


9 


1 , p, 1, 

/ X O I in —o— + — ill 

2 ( 47 r) 3 m B2 \ g 2 V2 2 2m 


+ 0(g 4 


?M2 ) 


(5.14) 


The condensate {4^4)^9yi2) can be measured non-perturbatively by subtracting proper 
counterterms [24] from lattice measurements extrapolated to the infinite-volume limit, and 
extrapolating subsequently to the continuum limit. A continuum extrapolation has only been 
carried out at an unphysically small Higgs mass [21] but cutoff effects have been seen to be 
modest, as long as we are not in the broken phase. Therefore we make use of lattice results at 
the physical Higgs mass [ 6 ] only for —0.1 < y < 0.2 in terms of the parameters in eq. (A.2). In 
order to extrapolate to higher or lower temperatures, perturbative expressions are employed; 
their values are discussed in appendix A. The procedure is illustrated in fig. 2. 


4 For gl i 2 corrections are known up to 2-loop order [23]. The other parameters of MSM read: 


- 


3<7b2 


16m 
2 


+ ■ 


2 2 
#E1$E2 


4 

#E1 


E2 4(m E1 + m E2 ) 16m E1 


W'3B — rn3B — J— (<?BimEl + 3gB2 ? UE2) + 


167T 

— 6ffEi5E2 ( -h l n 

V 4e m E1 + m E2 


8(4ti 


15 5 l 2 (^+ln 2m 


9 


+ ■ 


4 / 1 I , 9 

- 5ei — + In : 


3 

+ Io 


4e 


2m t 


(5.10) 


(5.11) 


These are needed if one wants to deduce eq. (5.7) directly from eq. (35) of ref. [18] by taking m E1 , m E2 m 3 . 
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Figure 2: The MS renormalized 3d condensate (cf. eq. (5.9)) according to lattice simulations [6] and 
perturbation theory (cf. appendix A), as a function of y from eq. (A.2). The lattice data corresponds 
to a fixed lattice spacing a, which induces errors particularly at low temperatures [21]. For the results 
of sec. 7 the phenomenological interpolation indicated with the thin dashed line is employed. 


5.4. Combined expression 

Summing together the contributions from eqs. (5.3), (5.6), (5.7) and (5.14), most of the 
1/e-divergences cancel, and we get 


2 u 2 Z m ((^ (/)) 


'J' 4 


~ ^ 2 (^ 1 A*) ~ 


4z/ 4 


(47r) 2 T 4 e 


+ 0(g 6 ), 


(5.15) 


where the finite part reads 

= — 7^r{l + 


v 


2 


1 - 


2(47t) 2 
3 


{g\ + 3g\ - 8A) ln(^) - Ah 2 In(^r) 


(^V)3d(5: 


2 ) 

M2 / 


3 T 2 1 2(4t r) 2 

^oTE 


(^+3 52 2 )(41n% + 31n i ^+ 7E + ^ +^1^ 




4h?ln(§^) + 8Aln(^) 
* V 8ttT ) V 4vrT ) 


+ 


2iA 


(47r) 3 T 2 
81 /' 


gfr u-Ei + 3g/m 

E2 _|_ gfT | 




/7 4 T / 35 2m 
+ ( XT + In 


16m B1 4(m E1 +m E2 ) 2m E2 \2A 


9m2 


8i/ 4 / /re 7E \ 

(4tt) 2 T 4 \AttTJ ' 


(5.16) 


Apart from the last term the result is //-independent up to 0(g 6 ). This /./-dependence as well 
as the divergence in eq. (5.15) cancel against terms from the vacuum subtraction, cf. sec. 6. 
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6. Vacuum subtraction and renormalization 


Suppose that we compute the bare vacuum pressure, p 0B , in a perturbative loop expansion: 
Pob = YIiPob- I n a gauge-fixed computation, the result is a function of the Higgs expectation 
value, which can likewise be determined order by order: v = . Even though v is 

gauge-dependent, p 0B is gauge-independent order by order in perturbation theory. Inserting 
v into p 0B , we can re-expand the result as 


Pob = Pm (v {0) + U (1) + ...) + Pob (V 0) + •••) + ••• = Pob (^ (0) ) + Pob 




+ ... 


( 6 . 1 ) 


where we made use of p^} = p^ B and [Por] / ( v ^) = 0 - The terms not shown are of 0{g 6 ) in 
our power counting. Therefore, it is sufficient for our purposes to compute p 0B up to 1-loop 
level and insert the tree-level Higgs vacuum expectation value v 2 = v 2 /A into the expression. 
Dropping the superscripts, we get 


Pob — 


^{v 2 + 8v 2 )v 2 - i(A + <5A)u 4 


+ 


3 m% 
32vr 2 

4 

. ™H 

647T 2 

The counterterms read 


- + ln^ 


5 

+ 6 


1 1 P 2 

—h In —73 

e mi 


+ 


+ 


3m 4 
647T 2 
3 m 4 
167T 2 


1 p 2 5 

7 + ln S5 + 6 
1 , / 3 
7 + ln ^ + 2 


v 2 =u 2 /\ 


+ 0(g 6 ) . (6.2) 


8v 2 = 


<5A = 


3z^ 2 


( 47 r) 2 e 


(4vr) 


9i + 3pl 


+ hi + 2A 


r s , i + ^gjgj + 3fl2 
16 


-&?(*?- 2A) + A(4A - «±M) 


whereas at the minimum the masses take the values 

(g 2 i+g 2 2 W 


g 2 v 2 


rn z = 


m 2 H = 2u 2 


hit 


w 4A ’ 4A 

Most divergences cancel, and the vacuum result becomes 

4p 0 „ __ 4 v 4 


rru = 


2A 


p'A 


— A 3 (T; p) + 


(4vr) 2 T 4 e ’ 


where 

A 3 (T;/2) = 


+ 


Av 4 


XT 4 (4vr) 2 T 4 

3 v 4 

+ 


In £ + 2 

,2 2 




2<?| 


' 4Ap 2 5 

l n 2 2 c 
g^f 6 


A 2 (16vr) 2 T 4 
3v 4 h 4 2A fi 2 3 
A 2 (4vr) 2 T 4 [ 11 h 2 ^ 2 + 2 


,2 \ 2 


+ ( 9i + g 2 ) 


In 


4Ap 2 


5 

+ - 


i.g 2 + g 2 )v 2 6 


+ 0{g 6 


(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 


(6.7) 
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In practice, it is not useful to evaluate eq. (6.7) directly, because our renormalization scale 
will be p ~ 7 tT, which would introduce large logarithms if inserted into eq. (6.7). Rather, all 
vacuum parameters are first evaluated at a scale p 0 = m z , and then evolved from p 0 to the 
thermal scale through RG equations. As an illustration, the 1-loop RG equations for the two 
dimensionful parameters appearing in our analysis read 


du 2 

t-dfi 

3i/ 2 / 

87T 2 \ 

9 2 i+ 39 2 2 h 2 2X 

4 

V 

o 

3 

i 



^ dp 

II 

OO I 

to| 



( 6 . 8 ) 

(6.9) 


from where we get A 3 (T;p) = 4 p on /T 4 . At the scale p = p 0 , the running couplings are 
expressed in terms of physical parameters through 1-loop relations specified e.g. in ref. [13]. 
(For many parameters a higher loop order could be extracted from more recent literature, 
but such corrections are smaller than thermal uncertainties, so for the sake of simplicity and 
reproducibility we make use of explicit 1-loop expressions.) 


7. Phenomenological results 

Summing together the terms from eqs. (4.6), (5.15), (6.6) as dictated by eq. (3.7), all 1/e- 
divergences cancel, and we obtain 

= A,(r) + A 2 (T;p) + A 3 (T;p) + 0(9«). (7.1) 

Because of a cancellation between A 2 and A 3 the result is formally independent of the 
renormalization scale p, even though in practice a residual /i-dependence is left over, as a 
reflection of unknown higher order corrections. We write p, = airT, and vary a in the range 
a £ (0.5...2.0) in order to get one impression on the corresponding uncertainty. 

The result of eq. (7.1) is accurate up to and including the order 0(g 5 ). It is well known 
from studies of the pressure of QCD, however, that certain odd orders show anomalously 
poor convergence. In particular, whereas the 0(g 2 ) correction to the pressure provides for a 
reasonable approximation, the 0(g 3 ) correction is far off. For A, the order 0(g 4 ) is related to 
the 0(g 2 ) correction to the pressure, and the order 0(g 5 ) to the 0(g 3 ) correction (cf. sec. 4). 
For this reason, the numerically most accurate estimate for A can probably be obtained by 
restricting to 0(g 4 ). A result corresponding to this accuracy is shown for the observable of 
eq. (7.1) in fig. 3. (For the pressure we display also the 0(g 5 ) result in fig. 4(left).) 

In order to obtain other thermodynamic functions, the boundary value needed for eq. (3.2) 
should be fixed. We do this on the low-temperature side, making use of the results of 
ref. [25] 5 and thereby setting p(T 0 )/Tg ~ 10.91 at T 0 = 100 GeV. On the high-temperature 

5 Tabulated results can be downloaded from www.laine.itp.unibe.ch/eos06/. 
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0( g 4 ) 



Figure 3: The trace “anomaly” from eq. (7.1). The grey band reflects variations of the renormalization 
scale in the range /2 = (0.5...2.0)7 tT. Contributions of 0(^ 4 /T 4 ) are seen to cancel at low temperatures. 


side, we expect to match to the results of ref. [18], after changing the overall renormalization 
condition to our eq. (2.6) and by making the changes listed in appendix B. Like for A, we 
make use of the result of 0(g 4 ), with the negative 0(g 5 ) QCD contribution expected to lead 
to an underestimate [26, 27]. (We have also experimented with an approximate 0(g 6 ) QCD 
contribution as estimated in ref. [25], finding a result which is in between the (D(g 4 ) and 
0(g 5 ) ones, confirming that the 0(g 5 ) result is most likely on the low side.) 

After choosing an initial condition, other thermodynamic functions are obtained as follows: 
the pressure p/T 4 from eq. (3.2); the energy density from e/T l = A + 3p/T, the entropy 
density s = p' from s/T 3 = A + 4 p/T 4 ] the heat capacity c = e' from c/T 3 = TA' + 7A + 
12 p/T 4 ; the equation-of-state parameter from w = p/e = 1/(3 + AT 4 /p); and the speed of 
sound squared from <? s = p'/e' = s/c. Some of these functions are conveniently parametrized 
through 


g e s{T) = 


e(T) 


s(T) 


c(T) 


7T 2 T 4 

’ etn ' 

~2t t 2 T 3 ' 

’ etU ' 

"2tt 2 T 3 " 

30 


45 


15 


(7.2) 


Results for all of these functions are shown in fig. 4. (The heat capacity and the speed of 
sound squared are the most “difficult” quantities, because they require taking a numerical 
temperature derivative from A.) 

It may be wondered why our new results do not match exactly the previous ones for high 
temperatures [18], even though both have been computed up to the same order in g. The main 
reason is that they involve different sets of higher-order corrections. In particular, eq. (7.1) has 


12 






















12.0 


11.0 


10.0 - 


9.0 


8 . 0 , 


0(g ) 


0(g ) 
0( g 5 l 


250 300 


77 GeV 



0.34 


0.33 


0(g ) 


0.32 - 


0.31, 


1/3 

- w 


250 300 


100 200 300 400 

77 GeV 


Figure 4: Left: The Standard Model pressure. Middle: g e g, h e g, as defined in eq. (7.2). 
The heat capacity (parametrized by z e ff) shows a narrow peak as is characteristic of a rapid 
crossover. Right: the equation-of-state parameter w and the speed of sound squared c 2 s . The 
grey bands reflect variations of the renormalization scale in the range £l = (0.5...2.0)7 tT. The 
low and high-T results correspond to ref. [25] (T<110 GeV) and ref. [18] (T> 250 GeV). 
The close-ups illustrate our estimates of theoretical uncertainties on the high-T side. 


been evaluated “strictly” to 0(g 4 ) or (D(g 5 ), because this has led to fairly simple expressions. 
In contrast, the high-temperature result is more cumbersome, including more terms (with 
non-zero corrections of 0(1), 0(g 2 ), 0(g 3 )) and complicated “soft” contributions, because 
the thermal mass parameters associated with the Higgs and gauge fields are of the same 
order. For evaluating such contributions the experience from QCD suggests that it is not 
worth re-expanding them in terms of the original couplings but rather to evaluate p M “as is”. 
We have followed separate procedures for the two regimes, because their difference permits 
for us to get another impression on the magnitude of unknown higher-order corrections. 

Among the various types of error estimates that we have made, that based on scale vari¬ 
ations is clearly a lower bound, because it only probes the magnitude of special types of 
corrections. The error estimate originating from the mismatches of various computations 
should in principle be a more reliable one, because all types of corrections are included; nev¬ 
ertheless it should still be treated as a lower bound. In practice, depending on the observable, 
one of the two gives a more conservative error estimate. In fig. 4 both estimates are shown; 
consequently we expect the theoretical uncertainty of our analysis to be on the percent level. 
Ultimately, the true accuracy can only be judged by carrying out a non-perturbative analysis 
of the observables that we have considered. 

With the help of the functions in eq. (7.2), the relationship of time and temperature in the 
Early Universe (assuming a flat geometry) can be expressed as 

3 /ITmpi dT _ y/ges(T)h eS (T) 

2 V vr 3 T 3 dt ~ ie ff (T) ' 1 j 
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It is seen that a peak in heat capacity (i.e. visible in fig. 4(middle), leads to a short 
period of slower temperature change [28], and correspondingly a mildly increased abundance 
of produced particles if a particle production process is under way, or a reduced density of 
weakly interacting relics disappearing through a co-annihilation process. 


8. Conclusions 

Drawing on existing multiloop computations [18] and lattice simulations within a dimension- 
ally reduced effective theory [6], and complementing these through new 1-loop, 2-loop and 
3-loop computations needed for determining the “trace anomaly” of the electroweak theory, 
we have estimated the basic thermodynamic functions that play a role for Standard Model 
thermodynamics at temperatures between 100 GeV and 300 GeV (cf. fig. 4). These results 
can be matched, for most observables with modest (~ 1%) discrepancies, to perturbative 
computations at low [25] or high [18] temperatures. 

One finding of our study is that despite the high temperatures considered, radiative cor¬ 
rections are larger than might naively be anticipated. Larger corrections also imply that 
uncertainties are less well under control. On the low-temperature side, we believe that the 
results of ref. [25] do contain uncertainties because the analysis of the electroweak sector 
was based on a low loop order and because QCD corrections, whose convergence is slow at 
finite temperatures, are substantial. On the high-temperature side, the QCD corrections 
continue to be an issue (cf. fig. 4(left)). As an example of the physical significance of these 
uncertainties, it may be noted from fig. 4 that the effective numbers of degrees of freedom 
remain below the canonical value 106.75 in the whole temperature range considered, and even 
decrease modestly as the temperature increases above 300 GeV. Even though there could be 
valid physics reasons for the decrease (such as that the effective Higgs mass parameter is very 
small across the crossover but then increases again), a similar tendency also originates from 
the 0(g A ) QCD contribution and is then an artifact of the truncation. (The 0(g 5 ) correction 
decreases the effective numbers of degrees of freedom further but turns the results into slowly 
increasing functions). Another peculiar feature of the O(g^) result is that the speed of sound 
squared c 2 is ~ 0.1% above 1/3 at T > 300 GeV, however this is again reversed by the 0(g 5 ) 
correction. In general, the estimated theoretical uncertainty of our results is on the percent 
level, and no conclusions should be drawn from features finer than this. 

In comparison with ref. [25], we find results for the effective numbers of degrees of freedom 
that are about 1% lower at T > 300 GeV. The reason is mostly due to negative electroweak 
radiative corrections proportional to g\ and h 2 , which were omitted in ref. [25]. 

Apart from radiative corrections, another important ingredient in our results is the 3d 
condensate {$<!>) 3d (cf. eq. (5.14)), measured non-perturbatively on the lattice and subse¬ 
quently converted to the MS scheme. It is interesting to note that apart from the known 
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strong correlation of (</>'</>) 3( j with the anomalous rate of baryon plus lepton number violation 
(cf. refs. [30, 6] and references therein), the same quantity also plays a role for other cos- 
mologically relevant observables, such as the production rate of non-relativistic right-handed 
neutrinos [31] or the relationship between lepton and baryon number densities [32], There¬ 
fore, it seems well motivated to improve on the existing measurements [6] by including the 
U(l) subgroup and by taking the continuum limit for a physical Higgs mass. In addition, 
it would be helpful to measure the susceptibility related to this condensate, so that taking 
a numerical temperature derivative, needed for estimating the heat capacity or the speed of 
sound squared, could be avoided. 

Our interpolations of all thermodynamic functions shown in fig. 4 can be downloaded, 
for a wide temperature interval, from www.laine.itp.unibe.ch/eosl5/. (In these results 
we switch from one regime to another with a temperature gap of 15-20 GeV in between.) 
Despite the remaining uncertainties, we hope that these results turn out to be helpful as a 
background equation of state in Dark Matter or Leptogenesis computations operating in this 
temperature range. 
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Appendix A. 3d condensate in perturbation theory 


We collect here results for the MS renormalized condensate (</>^) 3d defined within the MSM 
effective theory and evaluated at the scale jl = g (cf. eq. (5.9)). 

For 3> 0, a 3-loop perturbative expression for (</>^) 3d can be extracted from eq. (35) of 
ref. [18], after subtracting the contributions containing adjoint scalar fields (the latter need 
to be corrected as explained in appendix B). The expression reads 
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As illustrated in fig. 2 (dotted lines), the result agrees surprisingly well with lattice data as 
soon as y > 0.2. 

In the broken symmetry phase, i.e. for rhj <C 0, an explicit loopwise result is available 
up to 2-loop level, however it shows poor convergence [29]. Therefore we make use of a 
numerically determined value (referred to as the “Coleman-Weinberg (CW) method”) which 
includes a subset of higher-order corrections and has been tested against lattice simulations 
in refs. [29, 21]. The result is only available in the approximation g M1 = 0 but the corrections 
originating from g M1 are expected to be small. A comparison with lattice data is shown in 
fig. 2 (long-dashed line), with the low-temperature deviation expected to be reduced by a 
continuum extrapolation [21]. 


Appendix B. Differences with respect to ref. [18] 


Our results make use of expressions first derived in ref. [18], however there are a few technical 
points on which we disagree with this reference. As already discussed in the text, an important 
issue is that of the overall renormalization condition, eq. ( 2 . 6 ), which was imposed in a 
different form in ref. [18] (what was subtracted there was the pressure of a hypothetical 
zero-temperature “symmetric phase”). 

In addition, the following typographic errors have been detected in ref. [18]: in eq. ( 8 ), 
—3 gy —> —N c gy] in eq. (31), g 2 —> g 2 s /(An) 2 \ in eq. (33), g'^ 2 —> g 3 4 and 2 U 2 —> f h' 3 2 ; on the 
1st line of eq. (35), In In 2 (m 2 _^ 3 m 2 ) i/ 2 ; on the 5th row of eq. (45), §§ -> f§ In. 

As far as notation goes, we have replaced g 2 g%, g' 2 -» g 2 , g 2 - -> h 2 , A —>• p,, and 7 —>• 
Because of the trickiness of hypercharge assignments for a general representation, we choose 
directly Standard Model group theory factors: in the notation of refs. [18, 19], 

<i F —)• 2 , d A —> 3 , C F —>■ — , C A —> 2 , T f —> — , n s —> 1 , N c —> 3 . (B.l) 

Only the number of fermion generations, denoted by n G , is left free in our results. In order 
to simplify the expressions somewhat, we also substitute h 3 g%/A, h' 3 ^ g 2 / 4. 

We disagree with the square brackets in eq. (198) of ref. [18]. In the notation of ref. [18], 
this term arises from the three-dimensional integral 


T = 
new — 


1 


(B.2) 


Jpqr ( P 2 + ^D)[( r ~P ) 2 + m i](Q 2 + m l)[{ r - q ) 2 + w|]r 4 

Because of an infrared divergence, the integral needs to be carried out in the presence of 
dimensional regularization. The infrared divergence is powerlike so that, in the end, the limit 
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e —>• 0 can be taken. The square brackets in eq. (198) of ref. [18] correspond to the combination 
[...] = 24(47r) 3 777Q777 3 X new , and we believe that this combination has to be replaced as 
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Consequently, a round bracket on the third-but-last line in eq. (35) of ref. [18] needs to be 
changed as 
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